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ABSTRACT 

O Increasingly large populations of disk galaxies are now being observed at increas- 

^ ! ingly high redshifts, providing new constraints on our knowledge of how such galaxies 

evolve. Are these observations consistent with a cosmology in which structures form 
hierarchically? To probe this question, we employ SPH/N-body galaxy scale simula- 
tions of late-type galaxies. We examine the evolution of these simulated disk galaxies 
from redshift 1 to 0, looking at the mass-size and luminosity-size relations, and the 
^ ! thickness parameter, defined as the ratio of scale height to scale length. The structural 

parameters of our simulated disks settle down quickly, and after redshift z = 1 the 
galaxies evolve to become only slightly flatter. Our present day simulated galaxies are 
larger, more massive, less bright, and redder than at z = 1. The inside-out nature of 



OO 



the growth of our simulated galaxies reduces, and perhaps eliminates, expectations of 
q [ evolution in the size-mass relation. 

^! Subject headings: galaxies: evolution — galaxies: formation — galaxies: disks 

Q [ galaxies: structure — numerical methods 

CO 

•rn . 1. INTRODUCTION 

Galaxies are major structural components of our Universe, yet their formation and evolution 
remain an outstanding mystery of contemporary astrophysics. The latest generation of telescopes, 
including Keck, Subaru, the Very Large Telescope, Gemini, and the Hubble Space Telescope, pro- 
vide snapshots of galaxies at various redshifts. The properties of disk galaxies out to z = 1 have 
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been the subject of several recent studies. These properties include the Tully-Fisher relation (e.g. 
Vogt et al. 1996; Bizyaev & Mitronova 2002; Bohm et al. 2004), the luminosity-size, magnitude- 
size, and mass-size relations (e.g. Lilly et al. 1998; Simard et al. 1999; Bouwens & Silk 2002; 
Barden et al. 2005), the size distribution (e.g. Shen et al. 2003; Ravindranath et al. 2004; Fergu- 
son et al. 2004) and the disk "thickness" (Reshetnikov, Dettmar & Combes 2003, hereafter RDC, 
Elmegreen et al. 2005). 

The evolution of the luminosity-size and mass-size relations of disk galaxies remains a con- 
troversial issue. In particular, the interpretation of such evolution depends acutely on the selection 
biases of the surveys. Several studies have found an increase in the 5-band surface brightness of 
~ 1 mag arcsec^ 2 , out to redshift 1 (e.g. Schade et al. 1996; Roche et al. 1998; Lilly et al. 1998). 
Taking into account selection effects Simard et al. (1999) found minimal evolution with redshift, 
work that was supported by Ravindranath et al. (2004), who found that the luminosity-size relation 
evolves by less than 0.4 mag arcsec -1 between z — 1.25 and z — 0. Yet these groups invoke a 
population of high redshift, high surface brightness galaxies in their interpretation. Different in- 
terpretations of completeness have resulted in Trujillo & Aguerri (2004) and Barden et al. (2005) 
finding an evolution in the rest frame V-band surface brightness of ~ 0.8 and ~ 1.0 mag arcsec~ 2 
out to z ~ 0.7 and z ~ 1.0, respectively. Barden et al. (2005) also find that disks at a given absolute 
magnitude are bluer, with lower stellar mass to light ratios at z ~ 1, resulting in weak evolution in 
the relation between stellar mass and disk size. 

RDC selected 34 edge-on disk galaxies in the Hubble Deep Field with apparent diameters 
larger than 1 .3" and unperturbed morphology. They found an evolution between distant and local 
disk galaxies in their relative thickness or flatness, as characterized by the ratio of scale height h z 
to scale length hi. Their results indicate that disks at redshift z ~ 1 are smaller, in absolute value, 
than present day disk galaxies, and have a thickness ratio, h z /h\ larger by a factor ~ 1.5. support 
for such a thickening of disk galaxies at high redshifts also comes from Elmegreen et al. (2005), 
who inferred from their study of galaxies in the Hubble Ultra-Deep-Field that high redshift disk 
galaxies are thicker than local spirals by a factor of ~ 2. 

Studies such as these are clearly increasing our knowledge of the evolutionary stages of pop- 
ulations of disk galaxies, yet the path taken for galaxies to evolve between these stages remains 
unclear, and reliant on modeling. The general scenario of galaxy formation within the cosmo- 
logical context, as envisioned by White & Rees (1978), involves extended dark halos of galaxies 
forming hierarchically by gravitational clustering of non-dissipative dark matter. The luminous 
components form by a combination of the gravitational clustering and dissipative collapse. A gen- 
eral framework for the formation of disk galaxies was outlined in Fall & Efstathiou (1980). In 
their model, a uniformly rotating, homogeneous protogalactic cloud of gas begins to collapse as it 
decouples from the Hubble flow. This protogalaxy is endowed with angular momentum from tidal 
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torques driven by surrounding structure (Peebles 1969). It is assumed that the baryons destined 
to form the disk receive the same tidal torques as the dark matter before much dissipation occurs. 
If the collapse is smooth, then specific angular momentum is conserved and the gas forms a thin, 
rapidly rotating disk (e.g. Dalcanton, Spergel, & Summers 1997; Mo, Mao, & White 1998). 

Detailed direct modeling has highlighted that the above disk formation scenario is far from 
simple within a Cold Dark Matter (CDM) Universe . Complex baryonic processes involving mul- 
tiple gas phases and energy feedback, for example from quasars, supernovae, and stellar winds, 
have proven difficult to model, yet these processes appear to be crucial in the formation and evo- 
lution of disk galaxies within the context of hierarchical structure assembly. Early simulations 
employed thermal feedback, following Katz & Gunn (1991), augmented somewhat in Navarro & 
White (1993), who added a variable fraction of supernova energy as kinetic energy. This type of 
feedback has been popular up until recent times (e.g. Steinmetz & Navarro 2002; Abadi et al. 
2003), even though it was known to be highly inefficient, with thermal energy radiating away over 
timescales shorter than the dynamical timescales Katz (1992). The effects of this inefficient feed- 
back, in light of the the overcooling problem (White & Rees 1978), was the rapid formation of stars 
in the early collapsing halos. Subsequent mergers of these halos resulted in loss of angular mo- 
mentum to the dark halo, resulting in disk galaxies being dominated by stellar halo components, 
and significantly deficient in angular momentum. Weil, Eke, & Efstathiou (1998) showed that 
the suppression of early radiative cooling in numerical simulations can result in more realistically 
sized disks. The problem was tackled in Thacker & Couchman (2000) by turning off cooling for a 
fixed time in gas within the SPH smoothing kernel of a supernova event. Although lacking detailed 
insights into the multiphase gas processes involved in galaxy formation, these models succeeded in 
regulating star formation in small halos, and resulted in increased angular momentum in the simu- 
lated disk galaxy. In Brook et al. (2004a), we simulated galaxies using identical initial conditions 
with these two different feedback methods. The thermal feedback following Katz & Gunn (1991) 
resulted in galaxies which are dominated by a high mass, high metallicity spheroidal stellar compo- 
nent, with a low mass disk component. By regulating star formation in low mass building blocks, 
the adiabatic feedback model, which follows Thacker & Couchman (2000), formed galaxies with 
a dominant disk component, and a low mass, low metallicity stellar halo. Sommer-Larsen and 
collaborators have examined other implementations of feedback, and also shown the importance 
of regulating cooling at early epochs (Sommer-Larsen, Gelato, & Vedel 1999; Sommer-Larsen & 
Dolgov 2001; Sommer-Larsen, Gotz, & Portinari 2003). Meanwhile, Governato et al. (2004) high- 
lights the importance of resolution in simulating disk galaxies. Implementing these findings has 
allowed models such as Robertson et al. (2004), Okamoto et al. (2005) and Governato et al. (2005) 
to make progress in reproducing disk galaxies with properties approaching those of observed disks. 

Our chemo-dynamical models for simulating disk galaxies (which we will refer to as sGALS, 
to avoid confusion with real galaxies) have been successful in forming disk galaxies with a domi- 
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nant young, metal-rich disk component, with a metal poor, old stellar halo, and an old, intermediate 
metallicity thick disk component. The sGALS have many features observed in the Milky Way and 
local disk galaxies, including metallicities, abundances, ages, and color gradients of the stellar 
populations of the thin and thick disks, metallicities and abundances of the stellar halo, and the 
relation between luminosity and metallicity of the stellar halo (Brook et al. 2004a; Brook et al. 
2004b [hereafter BKGF]; Brook et al. 2005 [hereafter BGMK]; Renda, Gibson, & Kawata 2005). 
Here we report on the evolution of the disk component of such sGALS, specifically after redshift 
~ 1. We deliberately follow the evolution of our sGALS at times after which they have settled 
into a disk morphology, and trace their luminosity, scale length, scale height, and central surface 
brightness. We hope to provide insight into how the disk component of late-type galaxies evolve in 
their relatively "quiescent phase," during which mergers are only minor and star formation occurs 
almost exclusively in the disk regions. 

In this paper, we concentrate on the evolution of the size-mass and size-central surface bright- 
ness relations and the "thickness" of four simulated disk galaxies (sGALS), and what these tell 
us about interpreting observations. We also examine the evolution of the magnitude and color of 
our simulated galaxies. Our study will not allow direct comparison with the observational studies 
which compare distant galaxy populations to "like" local galaxy populations; it is not clear which 
galaxies are likely progenitors of today's galaxies. Yet, by providing evolutionary paths that disk 
galaxies can be expected to follow, our study will provide insights into the interpretation of such 
observations, in particular in how to interpolate between the two populations. 



2. THE CODE AND MODELS 

The four sGALS studied here are the same as were studied in BGMK. They are formed using 
our galactic chemo-dynamical evolution code, GCD+, which self-consistently models the effects 
of gravity, gas dynamics, radiative cooling, and star formation. Full details of GCD+ can be found 
in Kawata & Gibson (2003). All galaxies are formed with the Adiabatic Feedback Model from 
Brook et al. (2004a). In this model, gas within the SPH smoothing kernel of SNe II explosions is 
prevented from cooling. This adiabatic phase is assumed to last for the lifetime of the lowest mass 
star which ends as a SNe II, i.e. the lifetime of an 8-M star (100 Myr). In the AFM, the energy 
released by SNe la, which do not trace starburst regions, is fed back as thermal energy. 

We employ a semi-cosmological version of GCD+, where the initial conditions consist of 
an isolated sphere of dark matter and gas (Katz & Gunn 1991). This top-hat overdensity has 
an amplitude, Si, at initial redshift, Zi, which is approximately related to the collapse redshift, 
z c , by z c = 0.36<5j(l + zi) — 1 (e.g. Padmanabhan 1993). Small scale density fluctuations are 
superimposed on each sphere, parameterized by a 8 , seeding local collapse and subsequent star 
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formation. Longer wavelength fluctuations are incorporated by imparting a solid-body rotation 
corresponding to a spin parameter, A = J\E\^ 2 /GM^, to the initial sphere, where J is the total 
angular momentum, E the total energy, and M tot is the total mass of the sphere. This simplified 
model allows us to run a suite of simulations at high resolution, whilst retaining the most important 
features of full cosmological simulations. The choice of large values of A, and initial conditions in 
which no major merger occurs in late epochs (z < 1), ensures that disk sGALS are formed. We 
use an Einstein-deSitter CDM model; parameters relevant to all sGALS include f2 = 1, baryon 
fraction, £l b = 0.1, H = 50km s _1 Mpc -1 , star formation efficiency, c* = 0.05, and <r 8 = 0.5. 
This choice was made merely for convenience. As we showed in BGMK, the consequences of 
using this cosmology with the initial conditions of this study, rather than the now standard ACDM 
cosmology, are negligible. Essentially, prior to the redshift of collapse z c ~ 2, a ACDM universe 
is almost identical to a Einstein-de Sitter universe, while after that redshift, the system is so dense 
that the background cosmology becomes irrelevant. For a detailed proof, we refer the reader to 
BGMK. Other parameters varied slightly between the four models, but all were chosen to result in 
Milky -Way-type sGALS. Table 1, reproduced from BGMK, shows the values of total mass, A, and 
z c for the four simulations. Along with these variations in parameters, different random seeds are 
incorporated in the initial conditions of the sGALS, creating evolutionary diversity in our sample. 
We employed 3891 1 dark matter and 3891 1 gas/star particles for each model. 



3. RESULTS 

Our model results in the formation of late-type sGALS which have a dominant disk compo- 
nent, with a low mass, metal poor halo component (BKGF). Figure 1 shows the rest frame £?-band 
surface brightness profile of sGALl, both edge-on (upper panels) and face-on (lower panels), at 
three epochs: z = 0.7 (corresponding to a lookback time of 7.3 Gyrs) z = 0.5 (5.7 Gyrs), and 
z = 0. Figures 2-4 show corresponding plots for sGALS 2-4, shown at z = 1 (corresponding to 
lookback time of 8.5 Gyrs) z = 0.5 (z = 0.4 for sGAL2, lookback time 4.8 Gyrs), and z — 0. 
These epochs are chosen to be: (a) after the thick disk formation epoch, when the old thin disk 

Table 1 . Model Parameters for the 4 sGALS. 
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Zc 


M tot (M Q ) 
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sGALl 


1.8 


5xlO n 


0.0675 


sGAL2 


1.9 


5xlO n 


0.0675 


sGAL3 


2.0 


lxlO 12 


0.0600 


sGAL4 


2.2 


5xlO n 


0.0675 
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has begun to form; (b) at an intermediate time; and (c) the present, respectively. The different time 
chosen for the early epoch of sGALl (z = 0.7) was taken as that simulated galaxy is undergoing 
a significant merger at z = 1, and we need our sGALS to have a settled disk structure for a fair 
comparison with distant galaxies which have been selected for their disk features. Similarly, a 
satellite is being accreted to sGAL2 at the intermediate epoch, z ~0.5, so we chose a later inter- 
mediate time, z = 0.4 for this sGAL. The surface brightness is calculated by assuming that each 
star particle is a single stellar population, and using the SSPs of Kodama & Arimoto (1997), with a 
simple Salpeter initial mass function (IMF) parameterized by a power-law slope of 1.35 and lower 
and upper mass of 0.1 and 6OM , respectively. The photometric properties of sGALS discussed 
here are all intrinsic values, where we avoid discussing the effect of the dust absorption. 

Total £?-band magnitudes of the sGALS are found by integrating the luminosity densities 
displayed in these images. The results are shown in column 3 of Table 2. We see that the galaxies 
have an average M B ~ —21.1 at the earliest epoch (which we will refer to as z ~ 1), an average 
of —20.6 at the intermediate epoch (z ~ 0.5) and —19.7 at z — 0. The higher luminosity at earlier 
times is largely due to the enhanced star formation at the epoch during which the thick disk forms 
(see BKGF), the epoch which precedes the "early epoch" (z ~ 1) analyzed in this study. The rate 
of star formation is lower during the more quiescent later epochs, resulting in evolution toward 
fainter magnitudes at the present. 

We repeat the above process in both the rest frame V and I bands, and derive B — I and V — I 
colors. The B — I and V — I colors are shown in columns 5 and 6 of Table 2, respectively. The 
galaxies become redder during the relatively quiescent late epochs studied, from an average B — I 
(V - I) of 1.1 (0.7) at to 1.5 (0.9) at z ~ 0.5, and 1.7 (1.0) at z = 0. These redder colors are 
due to the effects of a decreasing star formation rate, aging stellar populations, and more metal rich 
stellar populations being born in the disk region. The redder colors are not sufficient to overcome 
the decrease in £?-band luminosity; the sGALS also become fainter in the redder bands. 

The £?-band surface brightness profiles for the four sGALS are plotted against radius in Fig- 
ure 5, at the three epochs shown in Figures 1-4. We use open stars for the early phase of evolution 
(z ~ 1), open squares for the intermediate phase (z ~ 0.5), and solid triangles for z = in this 
and all ensuing plots in order to allow evolution of the properties of the disk sGALS to be easily 
followed. Outside the central bulge regions, the disk components can be approximated reasonably 
by exponential fits. This allows us to derive disk scale lengths for each sGAL at each epoch. This 
is done in two stages. An approximate scale length h* is derived by fitting between 4 and 20 kpc. 
This allows us to then make fits between l.bh* (to exclude the bulge) and \h\. Numerous tests ver- 
ified that this technique gave robust results, which are shown column 9 in Table 2. Scale lengths 
increase from an average of 2.9 at z ~ 1 to an average of 3.5 at z ~ 0.5, and to 4.1 by z = 0. The 
implication is that the galaxies grow from the "inside-out." 



-7- 



We plot the surface brightness profiles of the four sGALS perpendicular to the disk in Figure 6. 
We exclude the bulge region when calculating how the surface brightness varies with distance from 
the disk plane. Again, the profiles are well fitted with exponentials, and we fit between and 2.5 
scale heights (essentially out to 1 .5 kpc in Figure 6) using similar techniques as employed to derive 
scale lengths. The derived scale heights are shown in column 8 of Table 2. Scale heights remain 
close to constant, with an average of 0.63 at z ~ 1 to an average of 0.62 at z ~ 0.5, and an 
average of 0.65 at z — 0. We derive the thickness of the disks, defined as the ratio of scale 
height to scale length, h z /h h and show the results column 10 in Table 2. On average, the galaxies 
are "thicker" at the earlier epochs, by a factor of 1.4 at z ~ 1. This is due to the no evolution 
of the scale hight and the increasing scale length with decreasing redshift. Our simulation does 
show little evolution of the scale height, because the thick disk forms first, followed by the thin 
disk formation, as discussed in BKGF If the thick disk comes from heating of the thin disk, we 
should expect an increasing scale height with decreasing redshift, and a slower evolution of the disk 
thickness. Further observations of the thickness evolution will help constrain thick disk formation 
scenarios. The evolution of the thickness in our sGALS (in the sense that the galaxies become 
thinner) is driven by the inside-out thin disk formation. At the intermediate epoch, the sGALS 
show a marginal increase in thickness, by a factor of 1.1 at z ~ 0.5; galaxies settle rapidly to disk 
morphologies after their last "significant" merger (i.e. a merger likely to alter the morphology of 
the galaxy). 

Numerical heating is an issue when considering vertical structure in numerical simulations, 
but we have reason to believe that our results are not affected significantly. Previous studies at 
similar resolution to ours have shown numerical vertical heating to be relatively low (Font et al. 
2001; Mayer 2004). Figure 2 of BGMK shows the vertical velocity dispersion versus age for 
four sGALS in this study, which demonstrates relatively little vertical heating during the disks' 
quiescent period; in the order of 20% over the last ~ 8 Gyrs. Even if this were totally attributed 
to numerical heating (unlikely as processes such as spiral structure and substructure heating will 
contribute), our results will not be greatly affected, because the vertical scale height is proportional 
to the square root of the velocity dispersion, i.e h z oc ^fa~ z (e.g. Bottema 1993). Our measured scale 
height at z = will be, in the worst case, overestimated by ~10%. As far as relevance to our study, 
the important point that we note is that any artificial heating will not affect our results qualitatively. 
In the absence of such effects, our conclusions that our simulated disk galaxies become flatter 
(h z /hi decreases) will hold. Our quantitative estimates would then be underestimating this effect. 
Having said this, higher resolution studies will increase the reliability of the quantitative estimates 

The way that the structural parameters, scale length and scale height, evolve relative to one 
another can be ascertained from Figure 7. Also plotted as crosses are observations of local disk 
galaxies from Schwarzkopf & Dettmar (2000). Plotted as open circles are observations of 34 edge- 
on disk galaxies at z~l, taken from RDC. The z ~ 1 disk galaxies of RDC are found to be smaller 
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and thicker than local galaxies, as determined by previous studies, such as Schwarzkopf & Dettmar 
(2000), de Grijs (1998), Kregel et al. (2002) and Bizyaev & Mitronova (2002). Our sGALS evolve 
in this plot primarily by increasing in scale length, with less evolution in scale height, resulting in 
galaxies evolving across this diagram, from left to right. 

Central surface brightness, again in the rest frame in the 5-band, is derived by extrapolating 
the exponential fits from Figure 5 to r = 0. Results are shown in column 7 of Table 2. We find 
that our sGALS are brighter at earlier epochs, with an average fi B = 21.0 at z ~ 1, compared 
with fi B = 21.8 at z ~ 0.5, and ji B = 22.7 at z — 0. This is related to the higher star formation 
rates at the earlier epochs in our sGALS, as shown in Figure 5 of BGMK. This is consistent with 
observations, which point to star formation rates being higher at earlier epochs (e.g. Madau et al. 
1996; Hammer et al. 1997; Flores et al. 1999; Juneau et al. 2005). We plot the 5-band central 
surface brightness against the flatness parameter, hi/h z (the inverse of the thickness parameter), in 
Figure 8. The sGALS have higher surface brightness and are thicker at earlier times, closer to the 
merging epochs when the thick disks formed. As the thin disk grows larger in the more quiescent 
later stages of their evolution, the sGALS become less bright, and also thinner. At each epoch, it is 
also apparent that brighter sGALS are thicker, although this trend is not as strong as the trend with 
time, because the sGALS we have studied are all relatively quiescent throughout the epoch during 
which evolution is followed in this study. This supports the conclusions of RDC, who plotted the 
relation between hi/h z and central surface brightness of z ~ 1 galaxies (their Fig .7), and found 
that the brightest objects have the thickest disk. They speculated that enhanced star formation and 
enhanced thickness are the results of ongoing interactions and mergings. 

We derive the stellar mass M* of each sGAL at each epoch using the friends-of-friends al- 
gorithm, and confirm that sensitivity to our choice of linking length is insignificant. Results are 
displayed in column 4 of Table 2. The inside-out nature of the growth of the thin disk in our 
sGALS is apparent in the relationship between the stellar mass and scale length, shown in Fig- 
ure 9. The mass and size of the sGALS grow together. Although we do not want to overemphasize 
the quantitative nature of our results, using just four relatively homogeneous L* sGALS, we note 
that the ratio M*/hi does not evolve, with the average at each epoch staying remarkably close to 
1O 1O M kpc -1 . Qualitatively, it is clear that the thin disks of our sGALS increase in size as they 
grow in mass, and the disk evolves along the mass scale length relation. 



4. SUMMARY AND CONCLUSION 

The evolution of individual galaxies cannot be observed directly, but lookback studies are 
providing increasing information on properties of populations of "like" galaxies, out to redshift 1 
and beyond. These studies remain difficult, in large part because it is unclear how to select which 
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Table 2. Derived parameters of our four sGALS at present (z = 0), intermediate (z ~ 0.5) and 
early (z ~ 1) epochs, in columns: (3) 5-band magnitude (4) Stellar Mass (5) B-I color (6) V-I 
color (7) £?-band central surface brightness (8) scale height (9) scale length (10) thickness 

parameter 



redshift 
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BCmag) 


M* (M ) 


B-I 


V-I 


(U_B(mag/arcsec 2 ) 


/if(kpc) 
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z = 


sGALl 


-19.5 


2.9xl0 10 


1.7 
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4.5 


0.16 


z = 


sGAL2 


-19.6 


3.0xl0 10 
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22.1 


0.59 


3.7 


0.16 


z = 


sGAL3 


-20.2 


6.6xl0 10 


1.8 


1.1 


21.4 


0.60 


3.5 


0.17 


z = 


sGAL4 


-19.6 
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0.15 


z = 


mean 


-19.7 
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0.65 
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0.16 


z = 0.5 


sGALl 
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galaxies should be compared at each epoch. In addition, there are large uncertainties involved in 
correcting for dust and cosmological dimming, especially in the case of edge-on disk galaxies. 

Several recent studies have compared the observed evolution of disk galaxy properties with 
semi-analytic and infall type models. For example, Barden et al. (2005) and Trujillo et al. (2005) 
show that one of the most basic assumptions of the semi-analytic approach, where disk size scales 
in proportion to the virial radius of the dark halo (Mo, Mao, & White 1998), is at odds with 
observations; the observed size evolution is weaker than the predicted R oc H(z)~ 2 / 3 , where H(z) 
is the Hubble parameter as a function of redshift. The infall model of Bouwens & Silk (2002) 
assumes a star formation rate at each radius and time, according to the local density (Schmidt 
1959), and the infall of metal-free gas. This model is able to match observations of the size 
evolution, at given luminosity, out to z ~ 1.7 (Trujillo et al. 2005). In order to explain the evolution 
of disk galaxies, the necessity of accounting for effects such as mass build up, star formation, 
energy feedback, and gas infall self-consistently within a CDM hierarchical framework, motivates 
us to confront observations with direct numerical simulations. 

We examine four simulated disk galaxies (sGALS), and do so in an idealized manner, free 
from the difficulties posed in observation. Our sGALS form thick disks at high redshift (z ~ 
1.5 — 2) during a period of mergers of gas rich "building blocks" (BKGF), and rapidly settle 
into a disk galaxy morphology by z ~ 1. We examine their subsequent evolution to the present, 
comparing various properties of the galaxies at high (~ 1), intermediate (~ 0.5), and zero redshift, 
where we find clear trends. Our main results are: 

• Scale heights of disk galaxies show no evolution after z ~ 1. Scale lengths increase slowly, 
resulting in flatter disks at later times. This is consistent with the observations of RDC, who 
find that disk galaxies at z ~ 1 are smaller than present-day galaxies, and have a thickness 
ratio, h z / hi larger by a factor of ~ 1.5. 

• The "inside-out" growth of disks is reflected in the scale length increasing roughly propor- 
tional to the stellar mass, resulting in little to no evolution of the mass-size relation. This 
is consistent with interpretations of Trujillo & Aguerri (2004); Barden et al. (2005); Trujillo 
et al. (2005). Previously, few direct numerical studies have probed these issues. In a study 
of sGALS of various morphologies in a fully self-consistent ACDM simulation, Sommer- 
Larsen, Gotz, & Portinari (2003) examined two disk galaxies, and found that while one 
clearly evolved inside-out, the other evolves outside-in. Samland & Gerhard (2003) examine 
the build up of disk galaxies through a spherical gas infall model, perhaps an idealization of 
the quiescent period of disk galaxy formation, and finds that disks form inside-out. 

• Disk sGALS become more massive, less bright, and redder between redshifts z ~ 1 and 0. 
The evolution of the luminosity of disk galaxies remains observationally controversial, with 
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selection effects largely determining interpretations (for a discussion see Barden et al. 2005). 
The mass build up in this epoch is relatively quiescent, with aging stellar populations and 
high metallicity stars born in the disk accounting for the galaxy reddening. 

• Rapid periods of star formation due to mergers and interactions can result in both thicken- 
ing of the disk, and galaxy brightening (see Figure 8). This is consistent with the results of 
RDC, as well as local observations which show that interacting galaxies are thicker (Reshet- 
nikov & Combes 1997), and that a correlation exists between central surface brightness and 
"thickness," in the sense that higher surface brightness galaxies are thicker (Bizyaev 2000; 
Bizyaev & Mitronova 2002; Bizyaev & Kajsin 2004). 

The quantitative results of our study are summarized in Table 2. These results are subject 
to uncertainties inherent in our models, such as star formation and feedback prescriptions. Also, 
our simulations do not cover the range of galaxy masses, sizes, morphologies, and histories of the 
large observed populations. We therefore feel that more modeling of a greater population of high- 
and low-redshift sGALS within a fully cosmological simulation is required in order to make more 
meaningful quantitative predictions. Having said this, clear trends have emerged in this study of 
the evolution of the properties of disk sGALS formed in a hierarchical, Cold Dark Matter Universe. 

We thank Agostino Renda for his help. Simulations were performed on Swinburne University 
Supercomputer, those of the Australian Partnership for Advanced Computing, and on "Purple- 
haze," the Supercomputer facility at Universite Laval. CB and HM thank the Canada Research 
Chair program and NSERC for support. BKG, DK & JB acknowledge the financial support of the 
Australian Research Council through its Discovery Project program. DK acknowledges the finan- 
cial support of the Japan Society for the Promotion of Science, through a Postdoctoral Fellowship 
for research abroad. 
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Fig. 1. — Rest frame £>-band luminosity plots for sGALl, plotted edge-on (upper panels) and 
face-on (lower panels) at redshifts z = 0.7 (left panels), z = 0.5 (middle panels), and z = (right 
panels). 
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Fig. 2. — Rest frame B-band luminosity plots for sGAL2, plotted edge-on (upper panels) and 
face-on (lower panels) at redshifts z — 1 (left panels), z = 0.4 (middle panels), and z = (right 
panels). 
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Fig. 3. — Rest frame £>-band luminosity plots for sGAL3, plotted edge-on (upper panels) and 
face-on (lower panels) at redshifts z — 1 (left panels), z = 0.5 (middle panels), and z = (right 
panels). 
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Fig. 4. — Rest frame B-band luminosity plots for sGAL4, plotted edge-on (upper panels) and 
face-on (lower panels) at redshifts z — 1 (left panels), z = 0.5 (middle panels), and z = (right 
panels). 



- 19- 



o 



1 1 1 1 1 1 1 1 


I I I I I I I I I I 

az=0.0 


4 


I I I I I I 


II 


11111 

az=0.0 


■ir 


□ z=0.5 








□ z=0.4 


□ 


*z=0.7 




o 

A* 




*z=1 


□ 






□ 






" A* 

□ 

▲ ▲ 


****** 


< 2 

m 

o 




□ 

*A § £ 
▲ [ 

▲ 1 


A 


** M A flfl 


sGAL 1 

I I I I I I I I 


1 1 1 1 1 1 1 1 1 




sGAL 2 

1 1 1 1 1 1 


1 1 1 1 


A A A 

I I I I 



5 10 15 20 

radius (kpc) 



5 10 15 20 

radius (kpc) 



o 



I I I I I I I I 


i i i i i i i i 
a Z =0.0 
□ z=0.5 




*z=1 










□ * 




- A □ ^ . 

A° 




* A , 

A D □ ~ 6 

* A 






A 






A A n *^v 

A y - 

^ S g 66 
*a2 




sGAL 3 

i i i i i i i i 


i i i i i i i I i 



5 10 15 20 

radius (kpc) 



O 



I I I I I I I I I I 


i i i i 
a Z =0.0 _ 


□ 


□ z=0.5 


A 


*z=1 






- □ 




A 








□ ^ 




- Afiftft 




□ 44 




ADD 








A 36 




_ AAA 


* A 


- sGAL 4 

i i i i i i i i i i 


l l l l 



5 10 15 20 

radius (kpc) 



Fig. 5. — Rest frame £?-band luminosity (log scale) versus radius for our 4 sGALS, plotted at z — 1 
(z = 0.7 for SGAL1, open stars), z = 0.5 (z = 0.4 for sGAL2, open squares) and redshifts z = 
(solid triangles). Outside the bulge regions, the profiles are fitted reasonably by exponentials. 



-20- 



2.5 



© A I- 

J 1.5 



60 
O 



0.5 



I I I I I I I I I I 


*z=0.0 




□ z=0.5 - 




*z=0.7 J 










: A a 






$ * : 

a - 


: sGAL 1 

i i i I i i i i 


i i i i i 



2.5 



0.5 1 1.5 

height (kpc) 



2.5 



J 1.5 



60 
O 



0.5 



I I I I 




i 


iz=0.0 








□ z=0.5 - 








ftZ=1 - 


ft 








: ▲ n 










□ 








▲ 




ft 






□ 


* : 






▲ 


□ 








▲ 








a : 


\ sGAL 3 

i i i I 




i 


i i i i 



0.5 1 1.5 

height (kpc) 



© j i- 
J 1.5 



60 
O 



0.5 



I I I I I I I 


III 


' lz= 


=b.b 






□ z= 


=0.4 I 


-* 




*z= 


=1 J 




e 








ft 


□ 
▲ 










□ 

▲ _ 






ft 




I sGAL 2 

I I I I I I I 


i 


I 


i i 



2.5 



J 1.5 



60 
O 



0.5 



0.5 1 1.5 

height (kpc) 



I I I I I I I I 


ii 

az= 


=b.b 




□ z= 


=0.5 - 




ftZ= 


=1 1 


▲ 






1 


1 








▲ 


j sGAL 4 

i i i i i i i i 


II 


i i 



0.5 1 1.5 

height (kpc) 



Fig. 6. — Rest frame £?-band luminosity (log scale) versus height above the disk plane for our four 
sGALS, plotted at the same epochs, using the same symbols, as Figure 5. The profiles are well 
fitted by exponentials. 
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Fig. 7. — Log (scale length) versus log (scale height) for our four sGALS, with symbols respond- 
ing to the three epochs as defined in Figures 1-4. Individual sGALS are joined by lines. Also 
plotted are observations of local disk galaxies from Schwarzkopf & Dettmar (2000) (crosses), and 
observations of edge-on disk galaxies at z ~ 1 from RDC (open circles). Our sGALS evolve 
primarily by increasing scale length, with little evolution in scale height; disk growth is inside-out. 
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Fig. 8. — The thickness of the sGALS, as defined by the ratio of scale height to scale length, 
plotted against the 5-band rest frame central surface brightness. We show the early (open stars), 
intermediate (open squares) and present (solid triangles) epochs, as defined in Figures 1-4. Indi- 
vidual sGALS are joined by lines. The sGALS have higher surface brightness and are thicker at 
earlier epochs. At each epoch, it also is apparent that higher surface brightness sGALS are thicker, 
although this trend is not as strong as the trend with time. 
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Fig. 9. — Log of the stellar mass of the sGALS against log scale length, with symbols representing 
the epochs as shown in previous diagrams. Individual sGALS are joined by lines. As the sGALS 
become more massive, they also become larger, growing from the inside-out. 



